Can Baryonic Features Produce the Observed 100ft 1 Mpc Clustering? 
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ABSTRACT 

We assess the possibility that baryonic acoustic oscillations in adiabatic models may 
explain the observations of excess power in large-scale structure on lOO/i -1 Mpc scales. 
The observed location restricts models to two extreme areas of parameter space. In either 
case, the baryon fraction must be large (Qb/^o ^ 0.3) to yield significant features. The 
first region requires Qq ^ 0.2h to match the location, implying large blue tilts (n ^ 1.4) 
to satisfy cluster abundance constraints. The power spectrum also continues to rise 
toward larger scales in these models. The second region requires Qq ~ 1, implying 
Qb well out of the range of big bang nucleosynthesis constraints; moreover, the peak 
is noticeably wider than the observations suggest. Testable features of both solutions 
are that they require moderate reionization and thereby generate potentially observable 
(~ 1/iK) large-angle polarization, as well as sub-arc-minute temperature fluctuations. 
In short, baryonic features in adiabatic models may explain the observed excess only if 
currently favored determinations of cosmological parameters are in substantial error or 
if present surveys do not represent a fair sample of lOO/i -1 Mpc structures. 
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1. Introduction 

As the study of large-scale structure has pushed to 
ever larger scales, several data samples have suggested 
the presence of excess power confined in a narrow re- 
gion around the 100/i _1 Mpc scale. The first such claim 



was the pencil-beam redshift survey of Broadhurst et al 



(1990), in which six concentrations of galaxies separated 
by a perio dic spaci ng of 1 28/t -1 Mpc w ere seen. Later 

~ 1992 



work (e.g. Bahcall 1991; Guzzo et al 



Willmcr 



t al. 1994) has confirmed that these overdensities are 



indeed part of extended structures rather than small- 
scale anomalies, and new pencil beams show similar be- 
havior (Broadhurst et al. 1995). More recently, the 2- 
dimensional power spectrum of the Las Campanas Red- 



shift Survey (Landy et al. 1996, hereafter LCRS) and 
the 3-dimensional power spectrum of rich Abell clus- 



ters (Einasto et al. 1997 and references therein) reveal 
a narrow peak at similar scales, k w 0.06/iMpc -1 and 
0.052/iMpc~ 1 respectively. 

Other data sets show anomalies on these scales, al- 
though they are unable to resolve a narrow feature. 
Three-dimensional reconstructions based on the angular 



correlations in the APM survey (Gaztanaga & Baugb 
1997) suggest a sharp drop in the power spectrum in 
this region. Finally, high-redshift C IV absorption 
lines in quasar spectra were found to be correlated on 
100/i _1 Mpc scales ( |Quashnock et al 1996); if this is due 
to large-scale structure, it indicates greater power than 
expected. Hence, several different lines of observational 
inquiry suggest excess power on 100/i -1 Mpc scales, per- 
haps in the form of a narrow peak at wavenumbers 
~0.05/i- 0.06/1 Mpc" 1 . 

Cosmological models based on collisionless dark mat- 
ter (e.g. cold dark matter), when combined with power- 
law initial power spectra, produce smooth power spec- 
tra at late times. Such models therefore cannot match 
the feature described above. However, if baryons are 
present in significant quantities, the coupling between 
them and the cosmic microwave background photons at 
redshifts z ;> 1000 produces acoustic oscillations near 
the 100/i -1 Mpc scale (see Eisenstein & Hu 1997, here- 
after EH97, and references therein). 

In this Letter, we consider whether acoustic features 
in adiabatic models can explain the narrow peak in the 
power spectrum at lOO/i^Mpc scales. This choice of 
location for a feature immediately restricts us to two 
rather extreme regions of parameter space. We examine 
each of these in turn, detailing their requirements and 
predictions. 



Throughout this paper, fl is the total density of 
matter relative to the critical density; Qb is that of the 
baryons. The power-law exponent of the initial power 
spectrum is denoted n (n = 1 is scale-invariant). The 
Hubble constant is written as 100/ikms _1 Mpc _1 . We 
assume a cosmological constant to make the universe 
flat; an open universe would have a power spectrum of 
identical shape but with a less favorable normalization. 

2. Constraints 

A cosmological model with cold dark matter and 
baryons exhibits a power spectrum with a broad global 
maximum (hereafter, the peak) at small wavenumbers 
k £ 0.05/i Mpc" 1 and a series of oscillations (hereafter, 
the bumps) at larger wavenumbers k > 0.05/iMpc -1 
(c.f. Fig. 4). Therefore, one may either attempt to as- 
sociate the peak or the first bump with the observed 
100/i _1 Mpc feature. This yields two disjoint areas of 
parameter space which we display in Figure 1 for two 
different values of h. We will now discuss these two 
regions separately. 

2.1. Low-Oq Region 

The region on the left in Figure 1 corresponds to 
placing the first bump in the region 0.045/iMpc _1 < 
k < 0.07/iMpc _1 . The bump shifts to smaller scales 
(higher k) as flo increases, as reflected in the left-right 
limits. We take the bump location to be the position of 
corresponding ma ximum in the oscillatory piece of the 
transfer function ( EH97 , eq. 25). 

The lower bound on the baryon fraction Qb / comes 
from the requirement that the amplitude of the b ump, 
as measured using the decomposition of EH97 , ex- 
ceeds 20% in power. Smaller oscillations would not 
explain the observations. The upper bound on the 
baryon fraction comes from requiring the bump am- 
plitude to be less than a factor of 1.6; presumably 
larger oscillations would have caused the second bump 
at k w 0.12/i Mpc" 1 to be detected ( |Peacock fc Doddg 
1994). Note that while lowering h from 0.8 to 0.5 causes 
the allowed region to shift unfavorably to even lower fig, 
increasing h to 1.0 only marginally relaxes the bound 
on f^o- 

Hence, one is restricted to a low value of Qq, approx- 
imately less than 0.2h. For h « 0.8, this does not dras- 



tically violate nucleosynthesis (e.g. Tytler et al. 1996). 
However, the moderate baryon fraction needed to pro- 
duce the bumps also causes a significant suppression of 
power at k ;> 0.02/iMpc" 1 . For a COBE-normalized 
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Fig. 1. — Parameter space available for a strong feature 
at fc « 0.05/iMpc . The f2o-f2b/fio plane is shown 
for h = 0.5 (dashed lines) and h = 0.8 (solid lines). 
Left region: a bump between 0.045/i and 0.07/iMpc -1 
with an amplitude between 20% and 160% as marked 
by filled dots at 20%, 40%, 80%, and 160%. Right re- 
gion: the peak at fc ;> 0.03/i Mpc^ 1 with an additional 
requirement on its prominence (positive half- width half- 
maximum <; 0.4 decades [c.f. Fig. ||). Nucleosynthesis 
constraints of ili,h 2 



0.024 (Tytler et al. 1996) are 



shown (BBNS). 



(Bunn fc White 1997) and scale- invariant initial spec- 



trum (n = 1), the resulting values of the fluctuations on 
the cluster scale <r$ are less than 0.5. This is far smaller 
than the value (;> 1.0 for these fio) needed to reproduce 
the abundance of galaxy clusters. 

Adding a significant blue tilt (n ;> 1.4) can increase 
<T8 enough to satisfy the cluster abundance constraint. 
We have taken the lowest value in the literature ( |Eke| 
et ^ 1996) to provide conservative lower bounds on 
n. We display this situation in Figure 2. In general, 
larger amplitude features must be balanced by larger 
tilts. Note that adding a tensor contribution to COBE 
or removing the cosmological constant will decrease the 
power spectrum normalization and in turn require even 



higher tilts ( [White fc Silkj 1996). 

Tilts of n ;> 1.3 are difficult to realize in infla- 
tionary models (Garcia-Bellido & Linde 1997). Em- 
pirical constraints depend entirely upon the range in 
wavenumber used to define the tilt. The limited range 



Fig. 2. — The value of as is shown as a function of 
baryon fraction for a COBE-normalized flat Slo = 0.15, 
h = 0.8 model. Several different values of tilt are shown 
(solid lines). The value of as required to match the 
present-day cluster abundance (dashed lines) is taken 
to be a 8 = 0.5£V°- 53+0 - 13n ° ( [Eke et al.| 1996) with 
25% variation to reflect errors here and in the COBE- 
normalization. The shaded region represents models 
that satisfy the cluster abundance and bump location 
while having J> 20% power enhancement. 



of scales available to COBE DMR allows only a weak 
constraint (k « 10~ 3 /i Mpc -1 , n <J 1.8; Gorski et al 



1996). Combining COBE with degree-scale CMB obser- 
vations (k w 10~ 2 h Mpc -1 ) limits the tilt more severely. 
However both constraints may be relaxed if the universe 
were reionized moderately early. Blue tilts extending to 
smaller scales are constrained by arcmi nutc-scal c CMB 
observations (fc ~ lh Mpc -1 , n <J 2; Vishniac 1987), 
the absence of spectral distortions from dissipation of 
acoustic waves af ter thermalization (fc ~ 10 4 Mpc~ 1 , 
n <^ 1.5; Hu et al. 1994), and limits on primordial black 

holes (fc ~ 10 15 Mpc" 1 , n £ 1.3; |Grecn fc Liddlc| 1997). 
In summary, strong blue tilts that extend from COBE 
to the smallest observable scales are ruled out, but be- 
tween COBE and cluster scales, the situation is less 
restrictive as the slope may decrease at smaller scales. 

The parameter space remaining to the low-ilo region 
after the peak location, peak amplitude, cluster abun- 
dance and tilt constraints are applied is shown as the 
shaded regions in Figs. and |. 
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Fig. 3. — The width of the peak of the power spectrum, 
in decades of wavenumber, for an Qq = 1, h = 0.5 
model of varying baryon fraction. Displayed are the full- 
width at half-maximum (FWHM) and the positive half- 
width (HWHM), defined as the range in k between the 
maximum of the spectrum and its half maximum in the 
direction of increasing k. We show n — 1 models (solid 
lines) and n = 1.4 models (dashed lines). Increasing 
h decreases the curves by a small amount. Requiring 
a HWHM < 0.4 (dotted lines) eliminates low baryon 
models. 

2.2. High-f2 Region 

The region on the right in Figure 1 corresponds to 
placing the peak at k > 0.03/iMpc~ . Although the 
lower limit is well below the observational region 0.05/i— 
0.06ft Mpc -1 , the peak in these models is sufficiently 
broad that the exact maximum need not lie directly on 
the preferred scale to yield an enhancement of power. 
Figure 1 assumes n = 1; adding a blue tilt shifts the 
peak to higher k, increasing the allowed region. The 
region for n — 1.4 and ft — 0.5 is very similar to that 
shown for ft = 0.8. 

The peak is generically much broader than the bump. 
As the baryon fraction increases, the high-fc side of the 
peak steepens significantly, giving rise to a prominent 
and asymmetric feature. Two statistics characterizing 
the width of the peak for an CIq = 1, h = 0.5 model 
are shown in Figure 3. Here one sees that the full- 
width at half-maximum (FWHM) always exceeds 0.85 
decades in k for n = 1 and 0.65 decades for n = 1.4. 



Similarly the range in k over which the power spectrum 
drops from its peak to its half-maximum in the high-fc 
direction (positive HWHM) always exceeds 0.3 decades 
for n = 1 and 0.25 decades for n — 1.4. Adding a 
blue tilt steepens the low-fc side of the peak, thereby 
decreasing the width. 

Hence, even for large baryon fractions, the peak 
may be too broad when compared with the narrow fea- 
ture observed in the cluster power spectrum, LCRS, or 
pencil-beams. However, the sharp break on the small- 
scale side of the maximum may be a sufficient departure 
from the usual low-baryon spectral shape as to allow 
these models to be statistically consistent with these 
observations. To reflect this situation, we place a lower 
limit on the baryon fraction in Figure [l] by requiring 
that the positive HWHM be less than 0.4 decades in k 
(c.f. Fig. @). 

3. Discussion 

In Figure 4, we show a representative example from 
each of the allowed regions and overlay them with ob- 
servational data sets. The top two curves show an 
n = 0.12, n b = 0.04, ft = 0.8, n = 1.6 model. 
The first bump is located near k = 0.06ft Mpc -1 and 
contains significantly more power than a zero-baryon, 
T = HqIi = 0.25 model (dashed line). The first bump 



is prominent and well-matched to the Einasto et al 
(1997) power spectrum; a similar model would fit the 
LCRS data. However, the peak at larger scales is yet 
higher, implying that power should continue to rise as 
we look toward larger scales. This is a generic feature of 
this region of parameter space — avoided only by enor- 
mous blue tilts (n J> 2.3) — and may well be incompati- 
ble with the turnover in the power spectrum suggested 
by the APM surv ey flBaugh fc Efstathiou| 1993; |Gaz 
tanaga fc Baugh 1997). Non-linearities would likely 



help to wash out the series of bumps at smaller scales 
(k > O.lftMpc" 1 ). 

The bottom two curves show an Oq = 1, ^fe = 0.4, 
h = 0.6, n = 0.95 model. Again, the model has ex- 
cess power on 100ft -1 Mpc scales relative to a T = 0.25 
model. Because of the high baryon fraction, this model 
does in fact produce the as needed to match the ilo = 1 
cluster abundance. Of course, the baryon density is 
in complete violation of bounds from nucleosynthesis 
(Tytler ct al. 1996). Due to its large width, the peak 
feature provides only a marginal, but perhaps adequate, 
fit to the Einasto et al. (1997) data. 



Although unusual, these models need not be at odds 
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Fig. 4. — Two representative models from the studied 
regions compared with large-scale structure data. The 
top two solid curves show an fio = 0.12, = 0.04, h = 
0.8, n — 1.6 model; the bottom two solid curves show 
an f^o = 1, Oft = 0.4, h = 0.6, n — 0.95 model. Dashed 
curves are a T = Q,oh = 0.25 zero-baryon model for 



comparison. Data sets are APM (Gaztanaga & Baugb 
1997) (crosses), the compilation of Peacock & Dodds 
(1 994) (open box es), and the cluster power spectrum 



of Einasto et al. (1997) (solid boxes). All error bars 
are 2-<r, and only data at k < 0.25/iMpc -1 have been 
plotted. All normalizations are arbitrary. The models, 
if COBE-normalized, partially reionized, and assumed 
flat, have as of 1.37 (no tensors) and 0.63 (with tensors), 
respectively. 

with current CMB observations. High baryon fractions 
tend to substantially enhance the first acoustic peak, 
and of course blue tilts enhance all power at smaller 
angular scales. If reionization were not invoked, the 
models would overproduce degree-scale anisotropics. 
However, with reionization corresponding to an opti- 
cal depth of r = 0.75 for the £lo = 0.12 model and 
t = 0.5 for the SIq = 1 model, the degree-scale predic- 
tions are suppressed to match current observations (e.g. 
Ncttcrfield et al. 1997). Because of the high baryon con- 
tent, these values of r correspond to rather low epochs 



Fig. 5. — CMB temperature anisotropy (including Vish- 
niac effect) and polarization predictions for the two 
models described in Figure 4. The models are consistent 
with the current observational limits (la error boxes on 
detections and 2a upper limits, see e.g. Smoot 1997 and 
references therein.) 



of reionization, z r ; = 33 and 13 respectively. Nor does 
the reionization overproduce secondary anisotropies; we 
find a Vishniac contribution ( Hu fc White 1996) across 
the ATCA band (£ w 4500) of AT/T = 2.7 x 10~ 6 and 
1.9 x 10 -6 respectively, well below the current limit of 
1.6 x 10~ 5 QSubrahmanyan et al| 1993). 



Two predictions of these models for the CMB (Fig. ||) 
are 1) that the second acoustic peak will be quite sup- 
pressed compared to the first and third for the high Q 
model, due to the high baryon content, and 2) that the 
high optical depth will produce substantial CMB po- 
larization levels, approaching bandpowers of 5 x 1Q — ' 
at COBE scales, as well as substantial sub-arc-minute 
temperature fluctuations. These are within reach of the 
current generation of CMB polarization experiments 
( Keating et al.| 1997) and interferometer experiments, 
respectively. 

In summary, adiabatic CDM+baryon universes with 
power-law initial power spectra produce the peak found 
at k w 0.05 — 0.06/iMpc _1 only in extreme regions of 
cosmological parameter space. Placing the first bary- 
onic bump at these wavenumbers requires values of £lo 
lower than those implied by dynamical mass measure- 
ments. This in turn requires extremely large blue tilts 



and moderate reionization. Avoiding tilts above n « 1.7 
necessitates a cosmological constant that exceeds limits 
from gravitational lens surveys QKochanck| l996). These 
models may also be in conflict with power spectrum ob- 



servations at even larger scales (k ~ 0.02/iMpc 
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On the other hand, placing the peak of the power 
spectrum at the observed scale requires high values of 
Qq. Such models need fib ;> 0.3 and even so provide a 
feature that is broader than the observations suggest. 
Dynamically favored values of Qqi sav ~0.3, place the 
first valley of the power spectrum at the desired place! 

The question remains as to whether the observations 
fairly sample the true power spectrum. The narrow 
width of the observed features may merely indicate that 
a small number of fc-modes are dominating the sam- 
ple. This is more likely if the distribution of ampli- 
tudes is non-Gaussian; for example, small nonlineari- 
ties in the density field increase the freque ncy of hot 
spots in realizations of the power spectrum (Amcndola 
1994). If the underlying theory has a broad peak around 
100/i -1 Mpc, different volumes may by chance produce 
spikes at slightly different locations, with the true width 
only being recovered in a larger survey. However, some 
underlying feature will still be required, as shown by 
the failure of simulations of trace-baryon models to re- 
produce the observations (LCRS). 

Can mildly non-linear evolution shift the location of 
the peak in the linear power spectrum? Second-order 
corrections to the real-space power spectrum (e.g. Jain 
& Bertschinger] 1994) act only to reduce the amplitude 
of features, although the effects are quite small on the 
scales in question. One possible loophole is coherent ef- 
fects in redshift space, which we plan to investigate us- 
ing the Zel'dovich approximation (Szalay et al. 1997). 
A second possibility is scale-dependent bias, for exam- 
ple if objects tend to trace the scale at which the power 
spectrum is steepest, rather than where it has its max- 
imum. 

Finally, one may consider models beyond those treated 
here. Isocurvature models (e.g. Peebles] 1987) produce 
a sequence of oscillations that are 90° out of phase 



with those of adiabatic models (Hu & Sugiyama 1996 



Sugiyama & Silk 1997). For CIq ~ 0.3, this places the 
peak of the power spectrum at the intended scale; the 
first bump is never relevant. Alternatively, one can 
place a feature directly in the initial power spectrum 
( Atrio-Barandela et al. 1997). Ongoing redshift surveys 
should measure the power spectrum to sufficient preci- 
sion to distinguish between these various explanations 
of the lOO/i -1 Mpc excess. 
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